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Chapter 1 
Anderson Transitions: Criticality, Symmetries, and Topologies 



Ch ! A. D. Mirlin*, F. Evers, I. V. Gornyit, and P. M. Ostrovsky* 

Institut fiir Nanotechnologie, 
Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany 



The physics of Anderson transitions between localized and metallic phases in 
disordered systems is reviewed. We focus on the character of criticality as well 
as on underlying symmetries and topologies that are crucial for understanding 
phase diagrams and the critical behavior. 



o . 

Q , 1.1. Introduction 

Quantum interference can completely suppress the difFusion of a particle in ran- 
^ ' doni potential, a phenomenon known as Anderson localization.^ For a given energy 

C*^ . and disorder strength the quantum states are either all localized or all delocalized. 

^ijT I This implies the existence of Anderson transitions between localized and metal- 

f~^ ■ lie phases in disordered electronic systems. A great progress in understanding of 

the corresponding physics was achieved in the seventies and the eighties, due to 
f^ ' the developments of scaling theory and field-theoretical approaches to localization, 

^^ . which demonstrated connections between the Anderson transition and conventional 

second-order phase transitions; see review articles^"^ and the book.'^ 

During the last 15 years considerable progress in the field has been made in 
K> , several research directions. This has strongly advanced the understanding of the 

^ ' physics of Anderson localization and associated quantum phase transitions and 

allows us to view it nowadays in a considerably broader and more general context.^ 
First, the symmetry classification of disordered systems has been completed. 
It has been understood that a complete set of random matrix theories includes, 
in addition to the three Wigner-Dyson classes, three chiral ensembles and four 
Bogoliubov-de Gennes ensembles.^ Zirnbauer has established a relation between 
random matrix theories, cr-models and Cartan's classification of symmetric spaces, 
which provides the mathematical basis for the statement of completeness of the 
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classification.* The additional ensembles are characterized by one of the additional 
symmetries - the chiral or the particle-hole one. The field theories (a-models) 
associated with these new symmetry classes have in fact been considered already in 
the eighties. However, it was only after their physical significance had been better 
understood that the new symmetry classes were studied systematically. 

Second, the classification of fixed points governing the localization transitions in 
disordered metals was found to be much richer than that of symmetries of random 
matrix ensembles (or field theories). The first prominent example of this was in 
fact given 25 years ago by Pruisken^ who showed that the quantum Hall transition 
is described by a a-model with an additional, topological, term. However, it is only 
recently that the variety of types of criticality — and, in particular, the impact of 
topology — was fully appreciated. Recent experimental discoveries of graphene and 
topological insulators have greatly boosted the research activity in this direction. 

Third, an important progress in understanding the statistics of wave functions 
at criticality has been made. Critical wave functions show very strong fluctuations 
and long-range correlations that are characterized by multifractality^'^'^^'^'^ imply- 
ing the presence of infinitely many relevant operators. The spectrum of multifractal 
exponents constitutes a crucially important characteristics of the Anderson transi- 
tion fixed point. The understanding of general properties of the statistics of critical 
wave functions and their mult ifr act ality was complemented by a detailed study - 
analytical and numerical - for a number of localization critical points, such as con- 
ventional Anderson transition in various dimensionalities, 2D Dirac fermions in a 
random vector potential, integer quantum Hall effect (IQHE), spin quantum Hall 
effect (SQHE), 2D symplectic-class Anderson transition, as well as the power-law 
random banded matrix (PRBM) model. 

Fourth, for several types of Anderson transitions, very detailed studies using 
both analytical and numerical tools have been performed. As a result, a fairly 
comprehensive quantitative understanding of the localization critical phenomena has 
been achieved. In particular, the PRBM model, which can be viewed as a ID system 
with long-range hopping, has been analytically solved on its critical line.^'^"'^'"'^^ This 
allowed a detailed study of the statistics of wave functions and energy levels at 
criticality. The PRBM model serves at present as a "toy model" for the Anderson 
criticality. This model possesses a truly marginal coupling, thus yielding a line 
of critical points and allowing to study the evolution of critical properties in the 
whole range from weak- to strong-coupling fixed points. Further recent advances 
in quantitative understanding of the critical behavior of Anderson transitions are 
related to exploration of network models of IQHE and its "relatives" from other 
symmetry classes, development of theories of disordered Dirac fermions, as well as 
large progress in numerical simulations. 

Finally, important advances have been achieved in understanding the impact of 
the electron- electron interaction on Anderson transitions. While this article mainly 
deals with non-interacting systems, we will discuss most prominent manifestations 
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of the interaction in Sec. 1.5 and 1.6.2. 

This article presents an overview of field with an emphasis on recent develop- 
ments. The main focus is put on conceptual issues related to phase diagrams, the 
nature of criticality, and the role of underlying symmetries and topologies. For a 
more detailed exposition of the physics of particular Anderson transition points and 
an extended bibliography the reader is referred to a recent review, Ref. 6. 

1.2. Anderson transitions in conventional symmetry classes 

1.2.1. Scaling theory, observables, and critical behavior 

When the energy or the disorder strength is varied, the system can undergo a 
transition from the metallic phase with delocalized eigenstates to the insulating 
phase, where eigenfunctions are exponentially localized,^ 

|V'(r)|~exp(-|r-ro|/e), (1.1) 

and ^ is the localization length. The character of this transition remained, however, 
unclear for roughly 20 years, until Wegner conjectured, developing earlier ideas of 
Thouless,^** a close connection between the Anderson transition and the scaling the- 
ory of critical phenomena. ^^ Three years later, Abrahams, Anderson, Licciardello, 
and Ramakrishnan formulated a scaling theory of localization,^® which describes 
the flow of the dimensionless conductance g with the system size L, 

d\ng/dlnL^ Pig). (1.2) 

This phenomenological theory was put on a solid basis after Wegner's discovery of 
the field-theoretical description of the localization problem in terms of a nonlinear 
cr-model,^'' Sec. 1.2.2. This paved the way for the resummation of singularities in 
perturbation theory at or near two dimensions^*'^^ and allowed to cast the scaling in 
the systematic form of a field-theoretical renormalization group (RG) . A microscopic 
derivation of the cr-model worked out in a number of papers^"^^^ has completed a 
case for it as the field theory of the Anderson localization. 

To analyze the transition, one starts from the Hamiltonian H consisting of the 
free part Hq and the disorder potential U{r): 

H = Ho + U{r) ; Ho = p^/2m. (1.3) 

The disorder is defined by the correlation function {U{r)U{r')); we can assume it 
to be of the white-noise type for definiteness, 

{U{r)U{r')) = {2TrpT)-^5{r - r'). (1.4) 

Here p is the density of states, r the mean free time and (...) denote the disor- 
der average. Models with finite-range and/or anisotropic disorder correlations are 
equivalent with respect to the long-time, long-distance behavior to the white noise 
model with rcnormalizcd parameters (tensor of diffusion coefficients).^'^ 
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The physical observables whose scahng at the transition point is of primary 
importance is the locahzation length ^ on the insulating side (say, E < Ec) and the 
DC conductivity a on the metallic side {E > Ec), 

^ o, [Ec - Ey , aoziE-E^y. (1.5) 

The corresponding critical indices v and s satisfy the scaling relation^^ s = v{d — 2). 
On a technical level, the transition manifests itself in a change of the behavior 
of the diffusion propagator, 

n(ri,r2;c.) - (Gi+^/2(ri,r2)G^_/2(r2,ri)), (1.6) 

where G^, G are retarded and advanced Green functions, 

G^'^(r,r') = (r|(ii;-i?±z77)-i|r'), r? ^ +0. (1.7) 

In the delocalized regime IT has the familiar diffusion form (in the momentum space) , 

Yi{ii,u)^2TTp{E)/{Dq^-iu:), (1.8) 

where p is the density of states (DOS) and D is the diffusion constant, related to 
the conductivity via the Einstein relation a — e^pD. In the insulating phase, the 
propagator ceases to have the Goldstone form (1.8) and becomes massive, 

n(ri,r2;c.)^^^(|ri-r2|/e), (1.9) 

with the function T{r) decaying exponentially on the scale of the localization length, 
J^{r/^) ^ exp(— r/^). It is worth emphasizing that the localization length ^ ob- 
tained from the averaged correlation function 11 = {G^G^), Eq. (1.6), is in gen- 
eral different from the one governing the exponential decay of the typical value 
Iltyp — exp(ln G^G"^). For example, in quasi-lD systems the two lengths differ by 
a factor of four.^^ However, this is usually not important for the definition of the 
critical index i'. We will return to observables related to critical fiuctuations of wave 
functions and discuss the corresponding family of critical exponents in Sec. 1.2.3. 

1.2.2. Field-theoretical description 

1.2.2.1. Effective field theory: Non-linear a -model 

In the original derivation of the CT-model,^^'^°"^^ the replica trick was used to per- 
form the disorder averaging. Within this approach, n copies of the system are 
considered, with fields (/)q,, a = l,...,n describing the particles, and the replica 
limit ri — > is taken in the end. The resulting c-model is defined on the n — > 
limit of either non-compact or compact symmetric space, depending on whether the 
fields 4>a are considered as bosonic or fermionic. As an example, for the unitary 
symmetry class (A), which corresponds to a system with broken time-reversal in- 
variance, the a-model target manifold is \J{n,n)/\J{n) x \J{n) in the first case and 
U(2n)/U(n) x U(n) in the second case, with n — > 0. A supersymmetric formulation 
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given by Efetov^ combines fermionic and bosonic degrees of freedom, with the field 
$ becoming a supervector. The resuhing cr-model is defined on a supersymmet- 
ric coset space, e.g. U(l, 1|2)/U(1|1) x U(l|l) for the unitary class. This manifold 
combines compact and non-compact features and represents a product of the hyper- 
boloid tf = U(l, 1)/U(1) X U(l) and the sphere S^ = U(2)/U(l) x U(l) "dressed" 
by anticommuting (Grassmannian) variables. While being equivalent to the replica 
version on the level of the perturbation theory (including its RG resummation), 
the supersymmetry formalism allows also for a non-perturbative treatment of the 
theory, which is particularly important for the analysis of the energy level and eigen- 
function statistics, properties of quasi-lD systems, topological effects, etc.^'^^'^"*'^^ 
Focusing on the unitary symmetry class, the expression for the propagator H, 
Eq. (1.6) is obtained as 

n{r^,r2:uj)^ J DQQ\l{r^)Ql\{r2)e-^^'^\ (1.10) 

where S[Q] is the cr- model action 

S[Q] ^If- I d'^r Str [-^(VQ)^ - 2iiokQ]. (1.11) 

Here Q — T~^AT is a 4 x 4 supermatrix that satisfies the condition Q^ = 1 and 
belongs to the u-model target space described above, A = diag{l, 1,-1,-1}, and 
Str denotes the supertrace. The size 4 of the matrix is due to (i) two types of the 
Green functions (advanced and retarded), and (ii) necessity to introduce bosonic 
and fermionic degrees of freedom to represent these Green's function in terms of a 
functional integral. The matrix Q consists thus of four 2x2 blocks according to 
its advanced-retarded structure, each of them being a supermatrix in the boson- 
fermion space. In particular, Qf^ is the boson-boson element of the RA block, and 
so on. One can also consider an average of the product of n retarded and n advanced 
Green functions, which will generate a cr-model defined on a larger manifold, with 
the base being a product of U(n,n)/U(n) x U(n) and U(2n)/U(n) x U(n) (these are 
the same structures as in the replica formalism, but now without the n — > limit). 
For other symmetry classes, the symmetry of the cr-model is different but the 
general picture is the same. For example, for the orthogonal class (AI) the 8x8 
Q-matrices span the manifold whose base is the product of the non-compact space 
0(2,2)/0(2) X 0(2) and the compact space Sp(4)/Sp(2) x Sp(2). The cr-model 
symmetric spaces for all the classes (Wigner-Dyson as well as unconventional) are 
listed in Sec. 1.3. 

1.2.2.2. RG in 2 + e dimensions; e-expansion 

The cr-model is the effective low-momentum, low-frequency theory of the problem, 
describing the dynamics of interacting soft modes - diffusons and cooperons. Its RG 
treatment yields a flow equation of the form (1.2), thus justifying the scaling theory 
of localization. The /3- function f3{t) = —dt/d\nL can be calculated perturbatively 
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in the coupling constant t inversely proportional to the diniensionless conductance, 
t — 1/2'KgJ^ This allows one to get the e-expansion for the critical exponents in 
2 + e dimensions, where the transition takes place at i* <C 1. In particular, for the 
orthogonal symmetry class (AI) one finds^ 

Pit) = et - 2*2 - 12C(3)t^ + 0{t^). (1.12) 

The transition point t, is given by the zero of the /3(i), 

i, ^e/2-(3/8)C(3)e^ + 0(e^). (1.13) 

The localization length exponent v is determined by the derivative 

V = -l/P'{U) = e-i - (9/4) C(3)e2 + 0{e^), (1.14) 

and the conductivity exponent s is 

s = i^e = l-(9/4)C(3)e3 + C'(e'^). (1.15) 

Numerical simulations of localization on fractals with dimensionality slightly above 
2 give the behavior of v that is in good agreement with Eq. (1.14).^* For the unitary 
symmetry class (A), the corresponding results read 

I3{t)^et-2t^ -Qt^ + 0{t'); (1.16) 

U = (e/2)^/2 - (3/2) (e/2)^/2 + 0{e^'^); (1.17) 

jy == l/2e - 3/4 + 0(e) ; s = 1/2 - (3/4)e + 0(6^). (1.18) 

In 2D (e — 0) the fixed point t* in both cases becomes zero: /3(i) is negative for any 
< > 0, implying that all states are localized. The situation is qualitatively different 
for the third — symplectic — Wigner-Dyson class. The corresponding /3-function is 
related to that for the orthogonal class via /3sp(0 — ^2/3o(— V^), yielding'' 

Pit) =€t + t^ - (3/4) C(3)t' + 0{t''). (1.19) 

In 2D the /3-function (1.19) is positive at sufficiently small t, implying the existence 
of a truly metallic phase at t < t*, with an Anderson transition at certain t* ^ 
1. This peculiarity of the symplectic class represents one of mechanisms of the 
emergence of criticality in 2D, see Sec. 1.4.1. The /3-functions of unconventional 
symmetry classes will be discussed in Sec. 1.3.5. 

1.2.3. Critical wave functions: Multifractality 

1.2.3.1. Scaling of inverse participation ratios and correlations at criticality 

Multifractality of wave functions, describing their strong fluctuations at criticality, 
is a striking feature of the Anderson transitions. ^^'■^^ Multifractality as a concept 
has been introduced by Mandelbrot.''^ Multifractal structures are characterized by 

'^For spinful systems, g here does not include summation over spin projections. 
''Here t = 1/ng, where g is the total conductance of the spinful system. 
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an infinite set of critical exponents describing tlie scaling of the moments of some 
distribution. This feature has been observed in various complex objects, such as 
the energy dissipating set in turbulence, strange attractors in chaotic dynamical 
systems, and the growth probability distribution in diffusion-limited aggregation. 
For the present problem, the underlying normalized measure is just |V'^(r)| and the 
corresponding moments are the inverse participation ratios (IPR) '^ 

P,= j d''Y\i,{v)\'">. (1.20) 

At criticality, Pq show an anomalous scaling with the system size L, 

(P,)-L'^(|V'(r)|2^)^L--«, (1.21) 

governed by a continuous set of exponents Tq. One often introduces fractal dimen- 
sions Dq via Tq—Dq{q — 1). In a metal Dq=d, in an insulator Dq~0, while at a 
critical point Dq is a non-trivial function of q, implying wave function multifractal- 
ity. Splitting off the normal part, one defines the anomalous dimensions A^, 

Tq = d{q - I) + Aq, (1.22) 

which distinguish the critical point from the metallic phase and determine the scale 
dependence of the wave function correlations. Among them, A2 = —rj plays the 
most prominent role, governing the spatial correlations of the "intensity" ji/'P, 

L2^(|V'(r)V'2(r')|)^(|r-r'|/i)-''. (1-23) 

Eq. (1.23) can be obtained from (1.21) by using that the wave function amplitudes 
become essentially uncorrelated at |r — r'| ^ L. Scaling behavior of higher order 
correlations, (|i/'^'^^(ri)V'^''^(r2) • • • '*/'^''"(rn)|), can be found in a similar way, e.g. 

id(9l+92) ^|^2gi (r^)^2q2 (^3)1) ^ L"^" -^'2 (jn - Ta |/i)^"+'2 '"^"^ '"^"^ . (1.24) 

Correlations of different (close in energy) eigenfunctions exhibit the same scaling, "^^ 
L^^{\^f{v)^^{r')\) \ f \r-r'\ Y^ 

where a; = e^ — ej, L^^ ^ {pto)^^/'^, p is the density of states, and |r — r'| < L^. 
For conventional classes, where the DOS is uncritical, the diffusion propagator (1.6) 
scales in the same way. 

In the field-theoretical language (Sec. 1.2.2), A^ are the leading anomalous di- 
mensions of the operators Tr((5A)'' (or, more generally, Tt:{QKY^ . . . Tr((5A)^™ with 
(ji+ . . . +qm=q)-'^'^ The strong multifractal fluctuations of wave functions at criti- 
cality are related to the fact that Ag<0 for (7>1, so that the corresponding operators 
increase under RG. In this formalism, the scaling of correlation functions [Eq. (1.23) 
and its generalizations] results from an operator product expansion. '^■^^'^''' 

'^Strictly speaking, Pq as defined by Eq. (1.20), diverges for sufficiently negative q {q < —1/2 for 
real «/> and q < —3/2 for complex ip), because of zeros of wave functions related to their oscillations 
on the scale of the wave length. To find Tq for such negative q, one should first smooth \^'^\ by 
averaging over some microscopic volume (block of several neighboring sites in the discrete version) . 
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1.2.3.2. Singularity spectrum f{a) 

The average IPR {Pq) are (up to the normahzation factor L'^) the moments of the 
distribution function 'P(IV'P) of the eigenfunction intensities. The behavior (1.21) 
of the moments corresponds to the intensity distribution function of the form 

n\r\) - j^^L-^^'^-"^^ (1-26) 

Indeed, calculating the moments (|i/'^'^|) with the distribution (1.26), one finds 

(P,) =L^(|i/^2g|^ _ f daL-i^+f'^''^ , (1.27) 



where we have introduced a — — In ji/i^l/lnL. Evaluation of the integral by the 
saddle-point method (justified at large L) reproduces Eq. (1.21), with the exponent 
Tq related to the singularity spectrum f{a) via the Legcndrc transformation, 

Tq ^ qa- f{a) , q = f'{a) , a ^ Tq. (1-28) 

The meaning of the function /(a) is as follows: it is the fractal dimension of the 
set of those points r where the eigenfunction intensity is |V'^(r)| ~ L~°'. In other 
words, in a lattice version of the model the number of such points scales as _L/(").36 

General properties of Tq and f{a) follow from their definitions and the wave 
function normalization: 

(i) Tq is a non-decreasing, convex function (r^ > 0, Tq < ), with tq — —d, 
Ti-0; 

(ii) f{a) is a convex function (/"(a) < 0) defined on the semiaxis a > with a 
maximum at some point ao (corresponding to q = under the Legendre transfor- 
mation) and /(ao) = d. Further, for the point ai (corresponding to g = 1) we have 
f{ai) = ai and f'{ai) = 1. 

If one formally defines /(a) for a metal, it will be concentrated in a single point 
a — d, with f{d) — d and f{a) — — oo otherwise. On the other hand, at criticality 
this "needle" broadens and the maximum shifts to a position ao > d, see Fig. 1.1. 

1.2.3.3. Symmetry of the multifractal spectra 

As was recently shown, ^^ the multifractal exponents for the Wigner-Dyson classes 
satisfy an exact symmetry relation 

A, =Ai_g, (1.29) 

connecting exponents with q < 1/2 (in particular, with negative q) to those with 
q > 1/2. In terms of the singularity spectrum, this implies 

f{2d-a) = f{a)+d-a. (1.30) 

The analytical derivation of Eqs. (1.29), (1.30) is based on the supersymmetric a- 
model; it has been confirmed by numerical simulations on the PRBM model^^'"^* 
(see Fig. 1.1b) and 2D Anderson transition of the symplectic class. '^^'^° 
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Fig. 1.1. Multifractality at Anderson transitions, a) Singularity spectrum f{a) in d = 2 + e 
for e = 0.01 and e = 0.2 (analytical), as well as 3D and 4D (numerical). With increasing d the 
spectrum gets broader, implying stronger multifractality. Inset: comparison between f{ct) for 3D 
and the one-loop result of the 2 + e expansion with e = 1 (solid). ^-"^ b) Singularity spectrum /(o) 
for the PRBM model. Evolution from weak to strong multifractality with decreasing parameter b 
is evident. Dashed lines represent /(2 — a) + a — 1, demonstrating the validity of Eq. (1.30).^* 

1.2.3.4. Dimensionality dependence of multifractality 

Let us analyze the evolution*^ from the weak-multifractality regime in d = 2 + e 
dimensions to the strong multifractality at d 3> 1. 

In 2 + e dimensions with e ^ 1 the multifractahty exponents can be obtained 
within the e-cxpansion, Sec. 1.2.2.2. The 4-loop results for the orthogonal and 
unitary symmetry classes read^^ 

A(0) = q{l - q)e + ^q(q - l)(q' - q + l)e^ + O(e^); 



Af) 



q{l-q)ie/2y/^--q^ 



; (1-31) 

{q-lfC{3)e^ + 0{e^/^). (1.32) 

and (1.32), we get the 



Keeping only the leading term on the r.h.s. of Eqs. (1.31) 
one-loop approximation for Tq which is of parabolic form. 

Numerical simulations^^ of the wave function statistics in 3D and 4D (Fig. 1.1a) 
have shown a full qualitative agreement with analytical predictions, both in the form 
of multifractal spectra and in the shape of the IPR distribution. Moreover, the one- 
loop result of the 2+e expansion with e—l describes the 3D singularity spectrum 
with a remarkable accuracy (though with detectable deviations). In particular, the 
position of the maximum, ao=4.03±0.05, is very close to its value ao=d+e implied 
by one-loop approximation. As expected, in 4D the deviations from parabolic shape 
are much more pronounced and q;o=6.5±0.2 differs noticeably from 6. 

The simulations*^ also show that fractal dimensions Dq = Tq/ {q— 1) with q^l 
decrease with increasing d. As an example, for q = 2 we have D2 ~ 2 — 2e in 2 -|- e 
dimensions, D2 = 1.3 ± 0.05 in 3D, and D2 — 0.9 ± 0.15 in 4D. This confirms the 
expectation based on the Bethe-lattice results (Sec. 1.2.5) that t^ ^ at rf ^ 00 
for q > 1/2. Such a behavior of the multifractal exponents is a manifestation of 
a very sparse character of critical eigenstates at d ^ 1 , formed by rare resonance 
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spikes. In combination with the relation (1.29) this imphes the hmiting form of the 
muhifractal spectrum at d — > oo, 

^ _/0' 9>l/2 .. „„. 

"«-\2d(,-l/2), q<l/2. ^'-^^^ 

This corresponds to f{a) of the form 

f{a) = a/2 , 0<a<2d , (1.34) 

dropping to — oo at the boundaries of the interval [0, 2d]. It was argued*-'^ that the 
way the mult ifr act ality spectrum approaches this limiting form with increasing d is 
analogous to the behavior found^"^ in the PRBM model with 6 <C 1. 

1.2.3.5. Surface vs. bulk multifractality 

Recently, the concept of wave function multifractality was extended*^ to the surface 
of a system at an Anderson transition. It was shown that fluctuations of critical 
wave functions at the surface are characterized by a new set of exponent r^ (or, 
equivalently, anomalous exponents A^) independent from their bulk counterparts, 

L'^-i(|^(r)p«)^L--.\ (1.35) 

T'g=d(q-l)+q^i+l + Al. (1.36) 

Here ^ is introduced for generality, in order to account for a possibility of non- 
trivial scaling of the average value, (|-0(r)P) oc L^"^^^, at the boundary in un- 
conventional symmetry classes. For the Wigner-Dyson classes, iJ, — 0. The nor- 
malization factor L'^^^ is chosen such that Eq. (1.35) yields the contribution of 
the surface to the IPR {Pq) = {J d'^r\tp{r)\'^''). The exponents A^ as defined in 
Eq. (1.36) vanish in a metal and govern statistical fluctuations of wave functions 
at the boundary, (|^(r)p'')/(|i/'(r)p)'' ^ L~ i, as well as their spatial correlations, 
e.g. L2(<i+M)(|^2(j.)^2(j.,)|^ _ (|r-r'|/L)^2. 

Wave function fluctuations are much stronger at the edge than in the bulk. As a 
result, surface exponents are important even if one performs a multifractal analysis 
for the whole sample, without separating it into "bulk" and "surface" , despite the 
fact that the weight of surface points is down by a factor 1/ L. 

The boundary multifractality was explicitly studied, analytically as well as nu- 
merically, for a variety of critical systems, including weak multifractality in 2D and 
2 + e dimensions, the 2D spin quantum Hall transition,''^ the Anderson transition 
in a 2D system with spin-orbit coupling,""^ and the PRBM model. ■^^ The notion of 
surface multifractality was further generalized'*^ to a corner of a critical system. 

1.2.4. Additional comments 

(i) For the lack of space we do not discuss the issues of IPR distributions at crit- 
icality and the role of ensemble averaging, as well as possible singularities in 
multifractal spectra, see the review 6. 
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(ii) Recently, an impressive progress was achieved in experimental studies of Ander- 
son transitions in various systems. '*^^'*^ The developed experimental techniques 
permit spatially resolved investigation of wave functions, thus paving the way to 
experimental study of multifractality. While the obtained multifractal spectra 
differ numerically from theoretical expectations (possibly because the systems 
were not exactly at criticality, or, in the case of electronic systems, pointing to 
importance of electron-electron interaction), the experimental advances seem 
very promising. 

(iii) Recent theoretical work^^ explains the properties of a superconductor-insulator 
transition observed in a class of disordered films in terms of multifractality of 
electronic wave functions. 

1.2.5. Anderson transition in d ^ oo; Bethe lattice 

The Bethe lattice (BL) is a tree-like lattice with a fixed coordination number. Since 
the number of sites at a distance r increases exponentially with r on the BL, it 
effectively corresponds to the limit of high dimensionality d. The BL models are 
the closest existing analogs of the mean-field theory for the case of the Anderson 
transition. The Anderson tight-binding model (lattice version of Eqs. (1.3), (1.4)) 
on the BL was studied for the first time in Ref. 50, where the existence of the 
metal-insulator transition was proven and the position of the mobility edge was de- 
termined. Later, the BL versions of the cr-model (1.11)^^'^^ and of the tight-binding 
model^*^ were studied within the supersymmetry formalism, which allowed to de- 
termine the critical behavior. It was found that the localization length diverges in 
the way usual for BL models, £, (x \E — Ec\~^ , where E is a microscopic parameter 
driving the transition. When reinterpreted within the effective-medium approxi- 
mation,"''^''''''' this yields the conventional mean-field value of the localization length 
exponent, ly ~ 1/2. On the other hand, the critical behavior of other observables is 
very peculiar. The inverse participation ratios Pq with q > 1/2 have a finite limit 
a,t E —>■ Ec when the critical point is approached from the localized phase and then 
jump to zero. By comparison with the scaling formula, Pq oc ^^^^ , this can be 
interpreted as Tg = for all q > 1/2. Further, in the delocalized phase the diffusion 
coefficient vanishes exponentially when the critical point is approached, 

D(y:n-^ln^n; n ^ exp{const\E - Ec\^^/^}, (1.37) 

which can be thought as corresponding to the infinite value, s — oo, of the critical 
index s. The distribution of the LDOS v = p{r)/{p) (normalized to its average 
value for convenience) was found to be of the form 

r{v)(y:n-^/^v-^^^ , n-^<$:v<$:n, (i.38) 

and exponentially small outside this range. Equation (1.38) implies for the moments 
of the LDOS: 

(w«) ex f^l«-i/2|-i/2, (139) 
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The physical reason for the unconventional critical behavior was unraveled in 
Ref. 56. It was shown that the exponential largeness of il reflects the spatial struc- 
ture of the BL: the "correlation volume" V^ (number of sites within a distance ^ 
from the given one) on such a lattice is exponentially large. On the other hand, 
for any finite dimensionality d the correlation volume has a power-law behavior, 
Vd{0 oc ^'^ oc |i5 — Ec\'"^, where ly ~ 1/2 at large d. Thus, the scale Q cannot ap- 
pear for finite d and, assuming some matching between the BL and large-d results, 
will be replaced by Vd((,)- Then Eq. (1.39) yields the following high-d behavior of the 
anomalous exponents A^ governing the scaling of the LDOS moments (Sec. 1.2.3), 

A,~d(l/2-|g-l/2|) , (1.40) 

or, equivalently, the results (1.33), (1.34) for the multifractal spectra r^, f{a). These 
formulas describe the strongest possible multifractality. 

The critical behavior of the conductivity, Eq. (1.37), is governed by the same 
exponentially large factor il. When it is replaced by the correlation volume Vd{^), 
the power-law behavior at finite d ^ 1 is recovered, a (x \E — £^c|^ with s ~ d/2. 
The result for the exponent s agrees (within its accuracy, i.e. to the leading order 
in d) with the scaling relation s ~ v{d — 2). 

1.3. Symmetries of disordered systems 

In this section we briefly review the symmetry classification of disordered systems 
based on the relation to the classical symmetric spaces, which was established in 
Refs. 7,8. 

1.3.1. Wigner-Dyson classes 

The random matrix theory (RMT) was introduced into physics by Wigner.^^ De- 
veloping Wigner's ideas, Dyson^^ put forward a classification scheme of ensembles 
of random Hamiltonians. This scheme takes into account the invariance of the 
system under time reversal and spin rotations, yielding three symmetry classes: 
unitary, orthogonal and symplectic. If the time-reversal invariance (T) is broken, 
the Hamiltonians are just arbitrary Hermitian matrices, 

H ^H'' , (1.41) 

with no further constraints. This set of matrices is invariant with respect to rota- 
tions by unitary matrices; hence the name "unitary ensemble" . In this situation, 
the presence or absence of spin rotation invariance (S) is not essential: if the spin 
is conserved, H is simply a spinless unitary-symmetry Hamiltonian times the unit 
matrix in the spin space. In the RMT one considers most frequently an ensemble 
of matrices with independent, Gaussian-distributed random entries - the Gaussian 
unitary ensemble (GUE). While disordered systems have much richer physics than 
the Gaussian ensembles, their symmetry classification is inherited from the RMT. 
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Let us now turn to the systems with preserved time-reversal invariance. The 
latter is represented by an antiunitary operator, T = KG , where C is the operator 
of complex conjugation and K is unitary. The time-reversal invariance thus implies 
H — KH^K^^ (we used the Hermiticity, H* = H^). Since acting twice with T 
should leave the physics unchanged, one infers that K*K — p, where p = ±1. As 
was shown by Wigner, the two cases correspond to systems with integer {p = +1) 
and half- integer {p = —1) angular momentum. If p = 1, a representation can be 
chosen where K ^ 1, so that 

H = H'^ . (1.42) 

The set of Hamiltonians thus spans the space of real symmetric matrices in this case. 
This is the orthogonal symmetry class; its representative is the Gaussian orthogonal 
ensemble (GOE). For disordered electronic systems this class is realized when spin 
is conserved, as the Hamiltonian then reduces to that for spinless particles (times 
unit matrix in the spin space). 

If T is preserved but S is broken, we have p = —I. In the standard repre- 
sentation, K is then realized by the second Pauli matrix, K = icTy, so that the 
Hamiltonian satisfies 

H = ayH'^Uy. (1.43) 

It is convenient to split the 2N x 2N Hamiltonian in 2 x 2 blocks (quaternions) in 
spin space. Each of them then is of the form q = qo(yo + i<li<7x+i<l2<7y + i<l30'z (where 
ao is the unit matrix and (Jx,y,z the Pauli matrices), with real g^, which defines a 
real quaternion. This set of Hamiltonians is invariant with respect to the group of 
unitary transformations conserving ctj,, U GyU'^ = ctj,, which is the symplectic group 
Sp(2Af). The corresponding symmetry class is thus called symplectic, and its RMT 
representative is the Gaussian symplectic ensemble (GSE). 

1.3.2. Relation to symmetric spaces 

Before discussing the relation to the families of symmetric spaces, we briefly remind 
the reader how the latter are constructed. ^^'^"^ Let G be one of the compact Lie 
groups SU(Af), SO{N), Sp(27V), and g the corresponding Lie algebra. Further, let 9 
be an involutive automorphism fl — > g such that 9^ = I but 9 is not identically equal 
to unity. It is clear that 9 splits g in two complementary subspaces, g — K(B^, 
such that 9{X) = X iov X € R and e{X) = -X for X e^. It is easy to see that 
the following Lie algebra multiplication relations holds: 

[i?,i?]ci?, [ii,V]c^, [*p,q2]cii. (1.44) 

This implies, in particular, that .^ is a subalgebra, whereas Cp is not. The coset space 
G/K (where K is the Lie group corresponding to ^) is then a compact symmetric 
space. The tangent space to G/K is *p. One can also construct an associated non- 
compact space. For this purpose, one first defines the Lie algebra g* ~ R Q) iCp, 
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which differs from g in that the elements in *p are multiphed by i. Going to the 
corresponding group and dividing K out, one gets a non-compact symmetric space 
G*/K. 

The groups G themselves are also symmetric spaces and can be viewed as coset 
spaces G x G/G. The corresponding non-compact space is G^ /G, where G^ is the 
complexification of G (which is obtained by taking the Lie algebra g, promoting it 
to the algebra over the field of complex numbers, and then exponentiating). 

The connection with symmetric spaces is now established in the following wayJ'^ 
Consider first the unitary symmetry class. Multiplying a Hamiltonian matrix by i, 
we get an antihermitean matrix X = iH. Such matrices form the Lie algebra u(A^). 
Exponentiating it, one gets the Lie group U(A^), which is the compact symmetric 
space of class A in Cartan's classification. For the orthogonal class, X = iH is 
purely imaginary and symmetric. The set of such matrices is a linear complement 
*P of the algebra M. = o{N) of imaginary antisymmetric matrices in the algebra 
Q = u{N) of antihermitean matrices. The corresponding symmetric space is G/K = 
\J{N)/0{N), which is termed AI in Cartan's classification. For the symplectic 
ensemble the same consideration leads to the symmetric space U(A^)/Sp(A^), which 
is the compact space of the class AIL If we don't multiply H hy i but instead 
proceed with H in the analogous way, we end up with associated non-compact 
spaces G* / K. To summarize, the linear space *P of Hamiltonians can be considered 
as a tangent space to the compact G/K and non-compact G* / K symmetric spaces 
of the appropriate symmetry class. 

This correspondence is summarized in Table 1.1, where the first three rows cor- 
respond to the Wigner-Dyson classes, the next three to the chiral classes (Sec. 1.3.3) 
and last four to the Bogoliubov-de Gennes classes (Sec. 1.3.4). The last two columns 
of the table specify the symmetry of the corresponding cr-model. In the supersym- 
metric formulation, the base of the cr-model target space A^b x Jvijr is the product 
of a non-compact symmetric space A^g corresponding to the bosonic sector and 
a compact ( "fermionic" ) symmetric space M.j^. (In the replica formulation, the 
space is A^g for bosonic or M.jr for fermionic replicas, supplemented with the limit 
n -^ 0.) The Cartan symbols for these symmetric spaces are given in the sixth col- 
umn, and the compact components Mjr are listed in the last column. It should be 
stressed that the symmetry classes of Alg and Mjr arc different from the symmetry 
class of the ensemble (i.e. of the Hamiltonian) and in most cases are also different 
from each other. Following the common convention, when we refer to a system as 
belonging to a particular class, we mean the symmetry class of the Hamiltonian. 

It is also worth emphasizing that the orthogonal groups appearing in the expres- 
sions for Aljr are O(A^) rather than SO(A^). This difference (which was irrelevant 
when we were discussing the symmetry of the Hamiltonians, as it does not affect 
the tangent space) is important here, since it infiuences topological properties of 
the manifold. As we will detail in Sec. 1.4,1.6, the topology of the a-model target 
space often affects the localization properties of the theory in a crucial way. 



Table 1.1. Symmetry classification of disordered systems. First column: symbol for the symmetry class of the Hamiltonian. Second 
column: names of the corresponding RMT. Third column: presence (+) or absence (— ) of the time-reversal (T) and spin-rotation (S) 
invariance. Fourth and fifth columns: families of the compact and non-compact symmetric spaces of the corresponding symmetry class. 
The Hamiltonians span the tangent space to these symmetric spaces. Sixth column: symmetry class of the c-model; the first symbol 
corresponds to the non-compact ("bosonic") and the second to the compact ( "fermionic" ) sector of the base of the cr-model manifold. 
The compact component Mp (which is particularly important for theories with non-trivial topological properties) is explicitly given in 
the last column. From Ref. 
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1.3.3. Chiral classes 

The Wigner-Dyson classes are the only allowed if one looks for a symmetry that is 
translationally invariant in energy, i.e. is not spoiled by adding a constant to the 
Hamiltonian. However, additional discrete symmetries may arise at some particular 
value of energy (which can be chosen to be zero without loss of generality) , leading 
to novel symmetry classes. As the vicinity of a special point in the energy space 
governs the physics in many cases (i.e. the band center in lattice models at half 
filling, or zero energy in gapless superconductors), these ensembles are of large 
interest. They can be subdivided into two groups - chiral and Bogoliubov - de 
Gennes ensembles - considered here and in Sec. 1.3.4, respectively. 

The chiral ensembles appeared in both contexts of particle physics and physics 
of disordered electronic systems about fifteen years ago.^^"^'"' The corresponding 
Hamiltonians have the form 



i.e. they possess the symmetry 



t,Ht, =. -H , (1.46) 



where Tz is the third Pauli matrix in a certain "isospin" space. In the condensed 
matter context, such ensembles arise, in particular, when one considers a tight- 
binding model on a bipartite lattice with randomness in hopping matrix elements 
only. In this case, H has the block structure (1.45) in the sublattice space. 

In addition to the chiral symmetry, a system may possess time reversal and/or 
spin-rotation invariance. In full analogy with the Wigner-Dyson classes, 1.3.1, one 
gets therefore three chiral classes (unitary, orthogonal, and symplectic). The cor- 
responding symmetric spaces, the Cartan notations for symmetry classes, and the 
cr- model manifolds are given in the rows 4-6 of the Table 1.1. 

1.3.4. Bogoliubov - de Gennes classes 

The Wigner-Dyson and chiral classes do not exhaust all possible symmetries of 
disordered electronic systems.^ The remaining four classes arise most naturally 
in superconducting systems. The quasiparticle dynamics in such systems can be 
described by the Bogoliubov-de Gennes Hamiltonian of the form 

N N 

H = J2 Kpcicp + - ^ (Aa044 - Kp^^o^cp) , (1.47) 

a/3 afi 

where c^ and c are fermionic creation and annihilation operators, and the N x N 
matrices h, A satisfy h ~ h^ and A-^ — —A, in view of hermiticity. Combining 
cjj, Ca in a spinor ifj^ — (cj^, Ca), one gets a matrix representation of the Hamiltonian, 
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H = tp'' Hip, where 

H=^_%_%), h = h^, A = -A^. (1.48) 

The minus signs in the definition of H result form the fermionic commutation 
relations between c^ and c. The Hamiltonian structure (1.48) corresponds to the 
condition 

H = -T^H^T^ (1.49) 

(in addition to the Hermiticity H = i?^), where t^ is the Pauli matrix in the particle- 
hole space. Alternatively, one can perform a unitary rotation of the basis, defining 
H = g^Hg with g = {1 + iTx)/V2. In this basis, the defining condition of class D 
becomes H = —H^, so that H is pure imaginary. The matrices X = iH thus form 
the Lie algebra so(2A^), corresponding to the Cartan class D. This symmetry class 
described disordered superconducting systems in the absence of other symmetries. 
Again, the symmetry class will be changed if the time reversal and/or spin 
rotation invariance are present. The difference with respect to the Wigner-Dyson 
and chiral classes is that now one gets four different classes rather than three. 
This is because the spin-rotation invariance has an impact even in the absence of 
time-reversal invariance, since it combines with the particle-hole symmetry in a 
non-trivial way. Indeed, if the spin is conserved, the Hamiltonian has the form 



N 



^/iy(c|^Cj^-Cj^4^) + Aij4^c]^^ + A*jCi^Cjt , (1.50) 



where h and A are NxN matrices satisfying h — h^ and A — A'^ . Similar to (1.48), 
we can introduce the spinors tpj = (cL,Ci^) and obtain the following matrix form 
of the Hamiltonian 

H=(t ^tV h = h\ A.= A^. (1.51) 



It exhibits a symmetry property 

H=-TyH^Ty. (1.52) 

The matrices H = iX now form the Lie algebra 5p(2N), which is the symmetry 
class C. 

If the time reversal invariance is present, one gets two more classes (CI and 
Dili). The symmetric spaces for the Hamiltonians and the cr- models corresponding 
to the Bogoliubov-de Gennes classes are given in the last four rows of the Table 1.1. 

The following comment is in order here. Strictly speaking, one should distinguish 
between the orthogonal group SO (TV) with even and odd N, which form different 
Cartan classes: S0(2A^) belongs to class D, while S0(2iV+l) to class B. In the 
conventional situation of a disordered superconductor, the matrix size is even due 
to the particlc-holc space doubling, see Sec. 1.3.4. It was found, however, that the 
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class B can arise in p-wave vortices/^ In the same sense, the class Dili should be 
split in Dlll-even and Dlll-odd; the last one represented by the symmetric space 
S0(4A^ + 2)/\J{2N + 1) can appear in vortices in the presence of time-reversal 
symmetry. 

1.3.5. Perturbative RG for cr-models of different symmetry classes 

Perturbative /3-functions for cr-models on all the types of symmetric spaces were 
in fact calculated^^'^^ long before the physical significance of the chiral and 
Bogoliubov-de Gennes classes has been fully appreciated. These results are im- 
portant for understanding the behavior of systems of different symmetry classes in 
2D. (We should emphasize once more, however, that this does not give a complete 
information about all possible types of criticality since the latter can be crucially 
affected by additional terms of topological character in the ct- model, see Sec. 1.4, 
1.6 below.) 

One finds that in the classes A, AI, C, CI the /3- function is negative in 2D 
in the replica limit (at least, for small t). This indicates that normally all states 
are localized in such systems in 2D. (This conclusion can in fact be changed in the 
presence of topological or Wess-Zumino terms. Sec. 1.4.1.) Above 2D, these systems 
undergo the Anderson transition that can be studied within the 2 + e expansion. 
Sec. 1.2.2.2. For the classes AIII, EDI, and CII (chiral unitary, orthogonal and 
symplectic classes, respectively) the /3(i) = in 2D, implying a line of fixed points. 
Finally, in the classes All, D, and Dili the /3-function is positive at small t, implying 
the existence of a metal-insulator transition at strong coupling in 2D. 

1.4. Criticality in 2D 

1.4.1. Mechanisms of criticality in 2D 

As was discussed in Sec. 1.2.2.2, conventional Anderson transitions in the orthogonal 
and unitary symmetry classes take place only if the dimensionality is d > 2, whereas 
in 2D all states are localized. It is, however, well understood by now that there is 
a rich variety of mechanisms that lead to emergence of criticality in 2D disordered 
systems. ^^ Such 2D critical points have been found to exist for 9 out of 10 symmetry 
classes, namely, in all classes except for the orthogonal class AI. A remarkable 
peculiarity of 2D critical points is that the critical conductance g* is at the same 
time the critical conductivity. We now list and briefly describe the mechanisms for 
the emergence of criticality. 

1.4.1.1. Broken spin-rotation invariance: Metallic phase 

We begin with the mechanism that has been already mentioned in Sec. 1.2.2.2 in 
the context of the Wigner-Dyson symplectic class (All) . In this case the /3-function 
[(1.19) with e = 0] is positive for not too large t (i.e. sufficiently large conductance), 
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SO that the system is metanic {t scales to zero under RG). On the other hand, for 
strong disorder (low t) the system is an insulator, as usual, i.e. /3{t) < 0. Thus, 
/3- function crosses zero at some i» , which is a point of the Anderson transition. 

This mechanism (positive /3-function and, thus, metallic phase at small t, with 
a transition at some i*) is also realized in two of Bogoliubov-de Gennes classes 

- D and DHL All these classes correspond to systems with broken spin-rotation 
invariance. The unconventional sign of the /3-function in these classes, indicating 
weak antilocalization (rather then localization), is physically related to destructive 
interference of time reversed paths for particles with spin s = 1/2. 

1.4.1.2. Chiral classes: Vanishing /3-function 

Another peculiarity of the perturbative /3-function takes place for three chiral classes 

- AIII, BDI, ad CII. Specifically, for these classes /3(i) = to all orders of the 
perturbation theory, as was first discovered by Gade and Wegner.^^'^"^ As a result, 
the conductance is not renormalized at all, serving as an exactly marginal coupling. 
There is thus a line of critical points for these models, labeled by the value of the 
conductance. In fact, the ct- models for these classes contain an additional term^^'^^ 
that does not affect the absence of renormalization of the conductance but is crucial 
for the analysis of the behavior of the DOS. 

1.4.1.3. Broken time-reversal invariance: Topological 9-term and quantum, 
Hall criticality 

For several classes, the a-model action allows for inclusion of a topological term, 
which is invisible to any order of the perturbation theory. This is the case when the 
second homotopy group 7r2 of the cr- model manifold M (a group of homotopy classes 
of maps of the sphere 5*^ into A^) is non-trivial. From this point of view, only the 
compact sector M.p (originating from the fermionic part of the supervector field) of 
the manifold base matters. There are five classes, for which t:2{M^f) is non-trivial, 
namely A, C, D, All, and CII. 

For the classes A, C, D the homotopy group n2{M.F) = Z. Therefore, the action 
S[Q] may include the (imaginary) 0-term, 

iStoplQ] = iSN[Q] , (1.53) 

where an integer N[Q] is the winding number of the field configuration Q{r). With- 
out loss of generality, 9 can be restricted to the interval [0, 2tt], since the theory is 
periodic in 9 with the period 2tt. 

The topological term (1.53) breaks the time reversal invariance, so it may only 
arise in the corresponding symmetry classes. The by far most famous case is the 
Wigner-Dyson unitary class (A) . As was first understood by Pruisken,^ the cr- model 
of this class with the topological term (1.53) describes the integer quantum Hall 
effect (IQHE), with the critical point of the plateau transition corresponding to 
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9=TT. More recently, it was understood that counterparts of the IQHE exist also 
in the Bogoliubov-de Gennes classes with broken time-reversal invariance - classes 
(^66-70 g^j-^j J) 71-74,79 fpjjgy -^vere called spin and t/ierma/ quantum Hall effects (SQHE 
and TQHE), respectively. 

1.4.1.4. Z2 topological term 

For two classes, All and CII, the second homotopy group is TT2{AiF) — ~^2- This 
allows for the 0-term but Q can only take the values and tt. It has been recently 
shown** that the cr-model of the Wigner-Dyson symplectic class (All) with a 0=7r 
topological angle arises from a model of Dirac fermions with random scalar potential, 
which describes, in particular, graphene with long-range disorder. Like in the case 
of quantum-Hall systems, this topological term inhibits localization. 

1.4.1.5. Wess-Zuvnino term 

Finally, one more mechanism of emergence of criticality is the Wess-Zumino (WZ) 
term that may appear in a-models of the classes AHI, CI, and DHL For these 
classes, the compact component M^p oi the manifold is the group H x HjH = H, 
where H is U(ri), Sp(2n), and 0(2n), respectively. The corresponding theories are 
called "principal chiral models" . The WZ term has the following form: 

^^wz(g) ^^y^Vy dse^,xStT{g-'d^g){g-'d,g){g-'dxg), (1.54) 

where k is an integer called the level of the WZW model. The definition (1.54) 
of the WZ term requires an extension of the a-model field g(r) = g{x,y) to the 
third dimension, < s < 1, such that g{r,0) = 1 and g(r, 1) — g{r). Such an 
extension is always possible, since the second homotopy group is trivial, tt2{H) = 0, 
for all the three classes. Further, the value of the WZ term does not depend on the 
particular way the extension to the third dimension is performed. (This becomes 
explicit when one calculates the variation of the WZ term: it is expressed in terms 
of g{r) only.) More precisely, there is the following topological ambiguity in the 
definition of S'wz(3). Since the third homotopy group is non-trivial, 1^^(11) — Z, 
Swzig) is defined up to an arbitrary additive integer n times 27rfc. This, however, 
does not affect any observables, since simply adds the phase nk x 27ri to the action. 
The WZ term arises when one bosonizes certain models of Dirac fermions^^ and 
is a manifestation of the chiral anomaly. In particular, a cr-modcl for a system of the 
AHI (chiral unitary) class with the WZ term describes Dirac fermions in a random 
vector potential. In this case the cr-model coupling constant is truly marginal (as 
is typical for chiral classes) and one finds a line of fixed points. On the other hand, 
for the class CI there is a single fixed point. The WZW models of these classes were 
encountered in the course of study of dirty d-wave superconductors'^'^* and, most 
recently, in the context of disordered graphene. We will discuss critical properties 
of these models in Sec. 1.4.2.3. 
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1.4.2. Disordered Dirac Hamiltonians and graphene 

Localization and criticality in models of 2D Dirac fcrmions subjected to various 
types of disorder have been studied in a large number of papers and in a variety 
of contexts, including the random bond Ising model, ®° the quantum Hall effect,*^ 
dirty superconductors with unconventional pairing, ^^"^^ and some lattice models 
with chiral symmetry. ^^ Recently, this class of problems has attracted a great deal 
of attention*^'^^"®^ in connection with its application to graphene. ^^'^^ 

One of the most prominent experimentally discovered features of graphene is 
the "minimal conductivity" at the neutrality (Dirac) point. Specifically, the con- 
ductivity^^"^^ of an undoped sample is close to e^//i per spin per valley, remaining 
almost constant in a very broad temperature range — from room temperature down 
to 30mK. This is in contrast with conventional 2D systems driven by Anderson 
localization into insulating state at low T and suggests that delocalization (and, 
possibly, quantum criticality) may emerge in a broad temperature range due to 
special character of disordered graphene Hamiltonian. 

In the presence of different types of randomness, Dirac Hamiltonians realize all 
ten symmetry classes of disordered systems; see Ref. 90 for a detailed symmetry 
classification. Furthermore, in many cases the Dirac character of fermions induces 
non-trivial topological properties (0-term or WZ term) of the corresponding field 
theory (cr-model). In Sec. 1.4.2.1 we review the classification of disorder in a two- 
flavor model of Dirac fermions describing the low-energy physics of graphene and 
types of criticality. The emergent critical theories will be discussed in Sec. 1.4.2.2- 
1.4.2.4. 

1.4.2.1. Symmetries of disorder and types of criticality. 

The presentation below largely follows Refs. 87,89. We concentrate on a two-flavor 
model, which is in particular relevant to the description of electronic properties 
of graphene. Graphene is a semimetal; its valence and conduction bands touch 
each other in two conical points K and K' of the Brillouin zone. In the vicinity 
of these points the electrons behave as massless relativistic (Dirac-like) particles. 
Therefore, the effective tight-binding low-energy Hamiltonian of clean graphene is 
a 4 X 4 matrix operating in the AB space of the two sublattices and in the K-K' 
space of the valleys: 

H = WoTgCrk. (1.55) 

Here T3 is the third Pauli matrix in the K~K' space, cr ~ {o'i,(T2} the two- 
dimensional vector of Pauli matrices in the AB space, and vq the velocity {vq ~ 10^ 
cm/s in graphene). It is worth emphasizing that the Dirac form of the Hamiltonian 
(1.55) does not rely on the tight-binding approximation but is protected by the 
symmetry of the honeycomb lattice which has two atoms in a unit cell. 

Let us analyze the symmetries of the clean Hamiltonian (1.55) in the AB and 
KK' spaces. First, there exists an SU(2) symmetry group in the space of the valleys. 
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Table 1.2. Disorder symmetries in graphene. The first five rows represent disorders pre- 
serving the time reversal symmetry To; the last four — violating Tq. First column: struc- 
ture of disorder in the sublattice {0-^) and valley (ri,) spaces. The remaining columns 
indicate which symmetries of the clean Hamiltonian are preserved by disorder.*^ 

structure Aj, Az Tp T± T^ Co Cx Cz CTg CT± CTz 

croTQ ++ + + + — — — — — — 

'^{1,2}T{1,2} -- +-- +-- +-- 

O-l,2T0 -+ +-+ +-+ +- + 

<^on,2 -- +-- --+ -- + 

CTSTS — + + — + — + — — + — 

o-3n,2 -- --+ --+ +-- 

0-OT3 -+ - + - - + - -I---I- 

0"1,2''"3 ++ — — — + + + — — — 

CT3T0 ++ --- --- + + + 



with the generators'^ 

Aa: = 0-3X1, Ay^(J3T2, A^=CroT3, (1-56) 

all of which commute with the Hamiltonian. Second, there are two more symmetries 
of the clean Hamiltonian, namely, time inversion operation (Tq) and chiral symmetry 
(Co). Combining Tq, Cq, and isospin rotations Aq, x.y.z, one can construct twelve 
symmetry operations, out of which four (denoted as T^) are of time-reversal type, 
four (C^) of chiral type, and four {CT^) of Bogoliubov-de Gennes type: 

To : A H^ (TiTiA (TiTi, Co '■ Ai-^ —a^ToAa^To, CTq : A h^ —a2TiA (T2T1, 

T^: A^ cr2ToA^cr2To, Cx- A^ -aQTiAaon, CT^ : A h^ -aiToA^aiTo, 

Ty : A^ <72T3A'^a2T3, Cy : A ^ -o-oT2A(ToT'2, CTy : A ^ -aiT^A^uiTz, 

Tj : A y~^ aiT2A^ aiT2, C^: A h^ -CT3T3ACT3T3, CT^ : A ^ -<T2T2A'^a2T2. 

It is worth recalling that the C and CT symmetries apply to the Dirac point {E ~ 0), 
i.e. to undoped graphene, and get broken by a non-zero energy E. We will assume 
the average isotropy of the disordered graphene, which implies that A^, and Ay 
symmetries of the Hamiltonian are present or absent simultaneously. They are thus 
combined into a single notation Aj^; the same applies to Tj_ and C±. In Table 1.2 
all possible matrix structures of disorder along with their symmetries are listed. 

If all types of disorder are present (i.e. no symmetries is preserved), the RG flow 
is towards the conventional localization fixed point (unitary Wigner-Dyson class 
A). If the only preserved symmetry is the time reversal (Tq), again the conventional 
localization (orthogonal Wigner-Dyson class AI) takes place.®'"' A non-trivial situa- 
tion occurs if either (i) one of the chiral symmetries is preserved or (ii) the valleys 
remain decoupled. In Table 1.3 we list situations when symmetry prevents local- 
ization and leads to criticality and non-zero conductivity at E=0 (in the case of 
decoupled nodes - also at nonzero E). Models with decoupled nodes are analyzed 
in Sec. 1.4.2.2, and models with a chiral symmetry in Sec. 1.4.2.3 (Co-chirality) and 
1.4.2.4 (C^-chirality). 
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Table 1.3. Possible types of disorder in graphene leading to criticality. The first three row cor- 
respond to Cz chiral symmetry leading to Gade-Wegner-type criticality, Sec. 1.4.2.4. The next 
three rows contain models with Co chiral symmetry (random gauge fields), inducing a WZ term 
in the tr-model action. Sec. 1.4.2.3. The last four rows correspond to the case of decoupled valleys 
(long-range disorder), see Sec. 1.4.2.2; In the last three cases the iT-modcl acquires a topological 
term with = n. Adapted from Ref. 89. 



Disorder 




Symmetries 


Class 


Criticality 


Conductivity 


Vacancies, strong potential impurities 
Vacancies + RMF 
o"3''"i,2 disorder 


Cz,To 

Cz 

Cz, Tz 


BDI 
AIII 

CII 


Cade 
Cade 
Cade 


^ Ae-'l-Kh 
sd 4eVvr/i 
Si ie^/TTh 


Dislocations 
Dislocations + RMF 




Co, To 
Co 


CI 
AIII 


WZW 
WZW 


ie-^/Txh 
ie'^/nh 


Ripples, RMF 




Az, Co 


2 X AIII 


WZW 


^e^l-Kh 


Charged impurities 
Random Dirac mass: 0-3 tq. 
Charged impurities -|- (RMF, 


0-3 r3 
ripples) 


Az, CTj_ 
Az 


2xAII 
2xD 

2xA 


9 = n 
9 = n 


40-5* or'' (4e^/7r/i)lnL 

4eV7r/i 

4<T£r 



''Numerical simulations®® reveal a flow towards the supermetal fixed point, u ~ {Ae'^ /nh) In L — > 00. 

1.4.2.2. Decoupled nodes: Disordered single-flavor Dirac fermions and 
quantum- Hall-type criticality 

If the disorder is of long-range character, the vaUey mixing is absent due to the lack 
of scattering with large momentum transfer. For each of the nodes, the system can 
then be described in terms of a single-flavor Dirac Hamiltonian, 



H = vo[cTk + crf,Vf,(r)]. 



(1.57) 



Here disorder includes random scalar (Vq) and vector (Vl.2) potentials and random 
mass (V3). The clean single-valley Hamiltonian (1.57) obeys the effective time- 
reversal invariance H = cr2H^a2- This symmetry (Ti) is not the physical time- 
reversal symmetry {Tq): the latter interchanges the nodes and is of no significance 
in the absence of inter-node scattering. 

Remarkably, single-flavor Dirac fermions are never in the conventional localized 
phase! More specifically, depending on which of the disorders are present, four 
different types of criticality take place: 

(i) The only disorder is the random vector potential (Vl,2)- This is a special case 
of the symmetry class AIII. This problem is exactly solvable. It is characterized by 
a line of fixed points, all showing conductivity Ae^/iih, see Sec. 1.4.2.3. 

(ii) Only random mass (V3) is present. The system belongs then to class D. 
The random-mass disorder is marginally irrelevant, and the system flows under RG 
towards the clean fixed point, with the conductivity Ae^ /irh. 

(iii) The only disorder is random scalar potential (Vq). The system is then in 
the Wigner-Dyson symplectic (AH) symmetry class. As was found in Ref. 88, the 
corresponding cr-model contains a Z2 topological term with 9 = i: which protects 
the system from localization. The absence of localization in this model has been 
confirmed in numerical simulations.^^ The scaling function has been found in Ref. 96 
to be strictly positive, implying a fiow towards the "supermetal" fixed point. 
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(iv) At least two types of randomness are present. All symmetries are broken in 
this case and the model belongs to the Wigner-Dyson unitary class A. It was argued 
in Rcf. 83 that it flows into the IQH transition fixed point. This is confirmed by the 
derivation of the corresponding cr-modelj*'*^'^^ which contains a topological term 
with = TT, i.e. is nothing but the Pruisken cr-model at criticality. A particular 
consequence of this is that the conductivity of graphene with this type of disorder 
is equal to the value ctJ} of the longitudinal conductivity a^x at the critical point of 
the IQH transition multiplied by four (because of spin and valleys) . 

If a uniform transverse magnetic field is applied, the topological angle becomes 
energy- dependent. However, at the Dirac point {E — 0), where a^y — 0, its value 
remains unchanged, 9 = tt. This implies the emergence of the half-integer quantum 
Hall effect, with a plateau transition point at E ^ 0. 



1.4.2.3. Preserved Cq chirality: Random gauge fields 

Let us consider a type of disorder which preserves the Co-chirality, H = —a^Ha^. 
This implies the disorder of the type cti. 2X9,1.2,3 being strictly off-diagonal in the a 
space. Depending on further symmetries, three different Co-chiral models arise: 

(i) The only disorder present is cri,2T3, which corresponds to the random abelian 
vector potential. In this case the nodes are decoupled, and the Hamiltonian decom- 
poses in two copies of a model of the class AIII. This model characterized by a line 
of fixed points has already been mentioned in Sec. 1.4. 2. 2. 

(ii) If the time-reversal symmetry Tq is preserved, only the disorder of the type 
'^i.2To.i,2 is allowed, and the problem is in the symmetry class CI. The model de- 
scribes then fermions coupled to a SU(2) non-abelian gauge field, and is a particular 
case of analogous SU(N) models. This theory flows now into an isolated fixed point, 
which is a WZW theory on the level k = -2Af.35.7r,97 

(iii) All Co-invariant disorder structures are present. This describes Dirac 
fermions coupled to both abelian U(l) and non-abelian SU(2) gauge fields. This 
model is in the AIII symmetry class. 

Remarkably, all these critical Co-chiral models are exactly solvable. In partic- 
ular, the critical conductivity can be calculated exactly and is independent on the 
disorder strength. A general proof of this statement based only on the gauge invari- 
ance is given in Ref. 87. (For particular cases it was earlier obtained in Refs.83,98). 
The critical conductivity is thus the same as in clean graphene. 



Ae^/irh. (1.58) 



Spectra of multifractal exponents and the critical index of the DOS can also be 
calculated exactly, see the review.^ 
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1.4.2.4. Disorders preserving C'z chirality: Gade-Wegner criticality 

Let us now turn to the disorder which preserves the C^-chirahty, H — —a^T^Ha^T^; 
according to Table 1.2, the corresponding disorder structure is cti. 2X0,3 and (70,3X1,2 If 
no time-reversal symmetries are preserved, the system belongs to the chiral unitary 
(AIII) class. The combination of Cz-chirality and the time reversal invariance Tq 
corresponds to the chiral orthogonal symmetry class BDI; this model has already 
been discussed in Sec. 1.4.1.2. Finally, the combination of Cz-chirality and T^- 
symmetry falls into the chiral symplectic symmetry class CII. The RG flow and 
DOS in these models have been analyzed in Ref. 81 In all the cases, the resulting 
theory is of the Gade-Wegner type.^^'^*^ These theories are characterized by lines 
of fixed points, with non- universal conductivity. It was found*''' "'^°° that for weak 
disorder the conductivity takes approximately the universal value, a ~ Ae'^/nh. In 
contrast to the case of Cq chirality, this result is, however, not exact. In particular, 
the leading correction to the clean conductivity is found in the second order in 
disorder strength.*^ 

1.5. Electron-electron-interaction effects 

Physically, the impact of interaction effects onto low-temperature transport and 
localization in disordered electronic systems can be subdivided into two distinct 
effects: (i) renormalization and (ii) dephasing. 

Renormalization. The renormalization effects, which are governed by virtual 
processes, become increasingly more pronounced with lowering temperature. The 
importance of such effects in diffusive low-dimensional systems was demonstrated 
by Altshulcr and Aronov, see Ref. 101. To resum the arising singular contributions, 
Finkelstein developed the RG approach based on the cr-model for an interacting 
system, see Ref. 102 for a review. This made possible an analysis of the critical 
behavior at the localization transition in 2+e dimensions in the situations when spin- 
rotation invariance is broken (by spin-orbit scattering, magnetic field, or magnetic 
impurities). However, in the case of preserved spin- rotation symmetry it was found 
that the strength of the interaction in spin-triplet channel scales to infinity at certain 
RG scale. This was interpreted as some kind of magnetic instability of the system; 
for a detailed exposition of proposed scenarios see Ref. 103. 

Recently, the problem has attracted a great deal of attention in connection with 
experiments on high-mobility low-density 2D electron structures (Si MOSFETs) 
giving an evidence in favor of a metal- insulator transition. ^'^'^ In Ref. 105 the RG 
for (j-model for interacting 2D electrons with a number of valleys A^ > 1 was an- 
alyzed on the two-loop level. It was shown that in the limit of large number of 
valleys N (in practice, A^ = 2 as in Si is already sufficient) the temperature of 
magnetic instability is suppressed down to unrealistically low temperatures, and a 
metal-insulator transition emerges. The existence of interaction-induced metallic 
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phase in 2D is due to the fact that, for a sufficiently strong interaction, its "delocal- 
izing" effect overcomes the disorder-induced locahzation. Recent works^^*'^^^ show 
that the RG theory describes weh the experimental data up to lowest accessible 
temperatures. We will see in Sec. 1.6 that the Coulomb interaction may also lead 
to dramatic effects in the context of topological insulators. 

The interaction-induced renormalization effects become extremely strong for cor- 
related ID systems (Luttinger liquids). While ID systems provide a paradigmatic 
example of strong Anderson localization, a sufficiently strong attractive interaction 
can lead to delocalization in such systems. An RG treatment of the correspond- 
ing localization transition in a disordered interacting ID systems was developed in 
Ref. 106, see also the book 107. Recently, the interplay between Anderson localiza- 
tion, Luttinger-liquid renormalization, and dephasing has been studied in detail in 
Ref. 120. 



Dephasing. We turn now to effects of dephasing governed by inelastic processes 
of electron-electron scattering at finite temperature T. The dephasing has been 
studied in great detail for metallic systems where it provides a cutoff for weak local- 
ization effects. ^°^ As to the Anderson transitions, they are quantum (zero-T) phase 
transitions, and dephasing contributes to their smearing at finite T. The dephasing- 
induced width of the transition scales as a power-law function of T. There is, how- 
ever, an interesting situation when dephasing processes can create a localization 
transition. We mean the systems where all states are localized in the absence of 
interaction, such as wires or 2D systems. At high temperatures, when the dephasing 
is strong, so that the dephasing rate tT (T) is larger than mean level spacing in 
the localization volume, the system is a good metal and its conductivity is given by 
the quasiclassical Drude conductivity with relatively small weak localization correc- 
tion.^"^ With lowering temperature the dephasing gets progressively less efficient, 
the localization effects proliferate, and eventually the system becomes an Anderson 
insulator. What is the nature of this state? A natural question is whether the 
interaction of an electron with other electrons will be sufficient to provide a kind 
of thermal bath that would assist the variable-range hopping transport, ^°® as it 
happens in the presence of a phonon bath. The answer to this question was given 
by Ref. 109, and it is negative. Fleishman and Anderson found that at low T the 
interaction of a "short-range class" (which includes a finite-range interaction in any 
dimensionality d and Coulomb interaction in d < 3) is not sufficient to delocal- 
ize otherwise localized electrons, so that the conductivity remains strictly zero. In 
combination with the Drude conductivity at high-T this implies the existence of 
transition at some temperature Tc. 

This conclusion was recently corroborated by an analysis^^"'^^^ in the framework 
of the idea of Anderson localization in Fock space. ^^^ In these works the temperature 
dependence of conductivity <t{T) in systems with localized states and weak electron- 
electron interaction was studied. It was found that with decreasing T the system 
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first shows a crossover from the weak-localization regime into that of "power-law 
hopping" over localized states (where ct is a power-law function of T), and then 
undergoes a localization transition. The transition is obtained both within a self- 
consistent Born approximation^ ^^ and an approximate mapping onto a model on the 
Bethe lattice. ^"'^° The latter yields also a critical behavior of a{T) above Tc, which 
has a characteristic for the Bethe lattice non-power-law form Incr(T) ^ (T — Tc)^'^ 
with K = 1/2, see Sec. 1.2.5. 

Up to now, this transition has not been observed in experiments'^, which indicate 
instead a smooth crossover from the metallic to the insulating phase with lowering 
y 113-116 rpj^g reason for this discrepancy remains unclear. An attempt to detect 
the transition in numerical simulations also did not give a clear confirmation of 
the theory, ^^^ possibly because of strong restrictions on the size of an interacting 
system that can be numerically diagonalized. On the other hand, a very recent 
work^^^ does report an evidence in favor of a transition of a Bethc-lattice character 
(though with different value of n) . 

1.6. Topological Insulators 

One of the most recent arenas where novel peculiar localization phenomena have 
been studied is physics of topological insulators. ^^^""'^^^ Topological insulators are 
bulk insulators with dclocalized (topologically protected) states on their surface.® 
As discussed above, the critical behavior of a system depends on the underlying 
topology. This is particularly relevant for topological insulators. 

The famous example of a topological insulator is a two-dimensional (2D) system 
on one of quantum Hall plateaus in the integer quantum Hall effect. Such a system 
is characterized by an integer (Chern number) n — ..., -—2, -—1, 0, 1, 2, ... which counts 
the edge states (here the sign determines the direction of chiral edge modes) . The 
integer quantum Hall edge is thus a topologically protected one-dimensional (ID) 
conductor realizing the group Z. 

Another (Z2) class of topological insulators^^^"^^^ can be realized in systems 
with strong spin-orbit interaction and without magnetic field (class AH) — and 
was discovered in 2D HgTe/HgCdTe structures in Ref. 125 (see also Ref. 127). A 
3D Z2 topological insulator^^^ has been found and investigated for the first time in 
Bii_a;Sba; crystals. Both in 2D and 3D, Z2 topological insulators are band insula- 
tors with the following properties: (i) time reversal invariance is preserved (unlike 
ordinary quantum Hall systems); (ii) there exists a topological invariant, which is 
similar to the Chern number in QHE; (iii) this invariant belongs to the group Z2 



"^Of course, in a real system, phonons are always present and provide a bath necessary to support 
the hopping conductivity at low T, so that there is no true transition. However, when the coupling 
to phonons is weak, this hopping conductivity will have a small prcfactor, yielding a "quasi- 
transition" . 

''Related topology-induced phenomena have been considered in Ref. 132 in the context of supcrfluid 
Hclium-3 films). 
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Table 1.4. Symmetry classes and "Periodic Table" of topological insula- 
tors. ^^'''^^^ The first column enumerates the symmetry classes of disordered 
systems which are defined as the symmetry classes Hp of the Hamiltonians (sec- 
ond column). The third column lists the symmetry classes of the classifying 
spaces (spaces of reduced Hamiltonians).^'^* The fourth column represents the 
symmetry classes of a compact sector of the sigma- model manifold. The fifth 
column displays the zeroth homotopy group 7ro(iJp) of the classifying space. 
The last four columns show the possibility of existence of Z and Z2 topological 
insulators in each symmetry class in dimensions d = 1,2,3,4. Adapted from 
Ref. 130. 

Symmetry classes Topological insulators 

p Hp Rp Sp TTojRp) d=l d=2 d=3 d=4 

Z Z 

Z2 Z 






AI 


BDI 


CII 


1 


BDI 


BD 


All 


2 


BD 


Dili 


Dili 


3 


Dili 


All 


BD 


4 


All 


CII 


BDI 


5 


CII 


C 


AI 


6 


C 


CI 


CI 


7 


CI 


AI 


c 



IL2 /L2 





Z2 Z2 Z 

Z 

Z 



'Li L2 



"-2 

Z 

z 



0' A AIII AIII 

1' AIII A A 



and reflects the presence or absence of delocalized edge modes (Kramers pairs). ^^^ 
Topological insulators exist in all ten symmetry classes in different dimensions, 
see Table 1.4. Very generally, the classification of topological insulators in d di- 
mensions can be constructed by studying the Anderson localization problem in a 
{d — l)-dimensional disordered system. ^'^^ Indeed, absence of localization of sur- 
face states due to the topological protection implies the topological character of the 
insulator. 

In Sec. 1.6.1 we overview the full classification of topological insulators and 
superconductors. ^^'''^^^ In Sec. 1.6.2 we discuss Z2 topological insulators belonging 
to the symplectic symmetry class All, characteristic to systems with strong spin- 
orbit interaction. Finally, in Sec. 1.6.3 we address, closely following Ref. 130, the 
interaction effects in Z2 topological insulators. 

1.6.1. Symmetry classification of topological insulators 

The full classification (periodic table) of topological insulators and superconduc- 
tors for all ten symmetry classes^'^ was developed in Refs. 135 and 134. This 
classification determines whether the Z or Z2 topological insulator is possible in 
the d-dimensional system of a given symmetry class. In this Section we overview 
the classification of topological insulators closely following Refs. 135 and 134, and 
discuss the connection between the classification schemes of these papers. 

All symmetry classes of disordered systems (see Section 1.3 and Table 1.1) can be 
divided into two groups: {A, AIII} and {all other}. The classes of the big group are 
labeled by p = 0, 1, . . . , 7. Each class is characterized by (i) Hamiltonian symmetry 
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class Hpi (ii) symmetry class Rp of the classifying space used by Kitaev;^'^^ (iii) 
symmetry class Sp of the compact sector Aip oi the sigma- model manifold. The 
symmetry class Rp of the classifying space of reduced Hamiltonians characterizes 
the space of matrices obtained from the Hamiltonian by keeping all eigenvectors 
and replacing all positive eigenvalues by +1 and all negative by —1. Note that 

Rp = -ffp+i, Sp — i?4-p. (1.59) 

Here and below cyclic definition of indices {0, 1, . . . , 7} (mod 8) and {0', 1'} (mod 
2) is assumed. 

For the classification of topological insulators it is important to know homotopy 
groups TTd for all symmetry classes. In Table 1.4 we list 7ro(i?p); other tt^ are given 

by 

nd{Rp)^TTo{Rp+d)- (1.60) 

The homotopy groups TTd have periodicity 8 (Bott periodicity). 

There are two ways to detect topological insulators: by inspecting the topology 
of (i) classifying space Rp or of (ii) the sigma- model space Sp. 

(i) Existence of topological insulator (TI) of class p in d dimensions is established 
by the homotopy group ttq for the classifying space Rp-d'- 

jTIofthetypeZ Jz 

JTIofthetypeZ. ^ -o(i?.-.) = |^^ (1-61) 

(ii) Alternatively, the existence of topological insulator of symmetry class p in d 
dimensions can be inferred from the homotopy groups of the sigma-model man- 
ifolds, as follows: 

JTIofthetypeZ ^^ ^d{Sp) = T 
\tI of the type Z2 ^^ ^d-i{Sp) = Z2 

The criterion (ii) is obtained if one requires existence of "non-localizable" bound- 
ary excitations. This may be guaranteed by either Wess-Zumino term in d — 1 
dimensions [which is equivalent to the Z topological term in d dimensions, i.e. 
TTd{Sp) = Z] for a QHE-type topological insulator, or by the Z2 topological term in 
d— 1 dimensions [i.e. itd-iiSp) = Z2] for a QSH-type topological insulator. 

The above criteria (i) and (ii) are equivalent, since 

T^diSp) = TTdiRi-p) = 7ro(i?4-p+<i)- (1.63) 

and 

, , fz forp = 0,4, , ^ 

MRp) = \^ ' ;; (1.64) 

[Z2 for p= 1,2. 

Below we focus on 2D systems of symplectic (All) symmetry class. One sees 
that this is the only symmetry class out of ten classes that supports the existence 
of Z2 topological insulators both in 2D and 3D. 
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1.6.2. Z2 topological insulators in 2D and 3D systems of class AH 



A Z2 class of topological insulators belonging to the symmetry class All was first 
realized in 2D HgTe/HgCdTe structures in Ref. 125. Such systems were found 
to possess two distinct insulating phases, both having a gap in the bulk electron 
spectrum but differing by edge properties. While the normal insulating phase has 
no edge states, the topologically nontrivial insulator is characterized by a pair of 
mutually time-reversed delocalized edge modes penetrating the bulk gap. Such state 
shows the quantum spin Hall (QSH) effect which was theoretically predicted in a 
model system of graphene with spin-orbit coupling. ^^^'^"^^ The transition between 
the two topologically nonequivalent phases (ordinary and QSH insulators) is driven 
by inverting the band gap.^^^ The Z2 topological order is robust with respect to 
disorder: since the time-reversal invariance forbids backscattering of edge states 
at the boundary of QSH insulators, these states are topologically protected from 
localization. 

For clean 2D QSH systems with a bulk gap generated by spin-orbit interaction, 
the Z2 invariant can be constructed from the Bloch wave functions on the Bril- 
louin zone^^^ and is somewhat similar to the Chcrn number in the standard QHE. 
Formally, if the Z2 index is odd/even there is an odd/even number m of Kramers 
pairs of gapless edge states (here m = is treated as even number) . In the pres- 
ence of disorder which generically back-scatters between different Kramers pairs, 
all the surface modes get localized if m was even in the clean system, while a single 
delocalized pair survives if m was odd. 

Disorder was found to induce a metallic phase separating the two (QSH and 
ordinary) insulators. ^^^'^^^ The transition between metal and any of the two insu- 
lators occurs at the critical value of conductivity 5 = 5* w 1.4; both transitions are 
believed to belong to the same universality class, see Sections 1.2.2.2 and 1.4.1.1. 
For g < g* all bulk states are eventually localized in the limit of large system, while 
for g > g* the weak antilocalization specific to the symplectic symmetry class drives 
the system to the "supcrmetallic" state, g — > 00. The schematic phase diagram for 
the nonintcracting case is shown in Fig. 1.2 (left panel). 

A related three-dimensional (3D) Z2 topological insulator was discovered in 
Ref. 126 where crystals of Bii_2;Sba; were investigated. The boundary in this 
case gives rise to a 2D topologically protected metal. Similarly to 2D topological 
insulators, the inversion of the 3D band gap induces an odd number of the surface 
2D modes. "'^^■^'"'^^^ These states in BiSb have been studied experimentally in Refs. 
126 and 128. Other examples of 3D topological insulators include BiTe and BiSe 
systems. ^^® The effective 2D surface Hamiltonian has a Rashba form and describes 
a single species of 2D massless Dirac particles (cf. Ref. 139). It is thus analogous to 
the Hamiltonian of graphene with just a single valley. In the absence of interaction, 
the conductivity of the disordered surface of a 3D topological insulator therefore 
scales to infinity with increasing the system size, see Section 1.4.2.2. 



July 7, 2010 0:7 World Scientific Review Volume - 9.75in x 6. Sin anderson-localization-SO. final 



Anderson Transitions 31 

1.6.3. Interaction effects on Z2 topological insulators of class All 

In this Section, we overview the effect of Coulomb interaction between electrons 
in topological insulators. ^^"^ Since a topological insulator is characterized by the 
presence of propagating surface modes, its robustness with respect to interactions 
means that interactions do not localize the boundary states. Indeed, arguments 
showing the stability of Z2 topological insulators with respect to interactions were 
given in Refs. 122,124,140 and 130. An additional argument in favor of persistence 
of topological protection in the presence of interaction is based on the replicated 
Matsubara sigma- model, in analogy with the ordinary QHE.^^^ This theory pos- 
sesses the same nontrivial topology as in the non-interacting case. 

Can the topologically protected 2D state be a supermetal (g ^ 00) as in the 
nonintcracting case? To answer this question the perturbative RG applicable for 
large conductivity g ^ 1 has been employed in Ref. 130. It is well known that 
in a 2D diffusive system the interaction leads to logarithmic corrections to the 
conductivity, ^°^ see Sec. 1.5. These corrections (together with the interference- 
induced ones) can be summed up with the use of RG technique. ^°^'^°'^ 

The one-loop equation for renormalization of the conductivity in the symplectic 
class with long-range Coulomb interaction and a single species of particles has the 
following form: 

/3(3) = ^ = -1/2. (1.65) 

Here —1/2 on the r.h.s. is a sum of the weak antilocalization correction 1/2 due 
to disorder and — 1 induced by the Coulomb interaction in the singlet channel. Ac- 
cording to Eq. (1.65), the negative interaction-induced term in /3(g) dominates the 
scaling at large g. Therefore, for g 3> 1 the conductance decreases upon renormal- 
ization and the supermetal fixed point becomes repulsive. 

Thus, on one hand, at 5 ^ 1 there is (i) scaling towards smaller g on the 
side of large g. On the other hand, surface states are topologically protected from 
localization, which yields (ii) scaling towards higher g on the side of small g. The 
combination of (i) and (ii) leads unavoidably to the conclusion that the system 
should scale to a critical state {g ~ 1). Indeed, there is no other way to continuously 
interpolate between negative (i) and positive (ii) beta functions: at some point I3{g) 
should cross zero. As a result, a critical point emerges due to the combined effect of 
interaction and topology. ^■^'^ In other words, if the system can flow neither towards 
a supermetal (g — > 00) nor to an insulator {g — )■ 0) it must flow to an intermediate 
fixed point {g ^ 1). Remarkably, the critical state emerges on the surface of a 3D 
topological insulator without any adjustable parameters. This phenomenon can be 
thus called "self-organized quantum criticality".^'^^ 

Let us now return to 2D Z2 topological insulators. The 2D disordered QSH 
system contains only a single flavor of particles, A^ = 1. Indeed, the spin-orbit 
interaction breaks the spin-rotational symmetry, whereas the valleys are mixed by 
disorder. As a result, the supermetal phase does not survive in the presence of 
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(a) no interaction 



(b) with interaction 
critical 




inverted normal 

band gap 



nverted normal 

band gap 



Fig. 1.2. The phase diagrams of a disordered 2D system demonstrating the QSH effect. Left: 
noninteracting case. Right: interacting case (with Coulomb interaction not screened by external 
gates). Interaction "kills" the supermetallic phase. As a result, the two insulating phases are 
separated by the critical line. Adapted from Ref. 130. 

Coulomb interaction: at g 3> 1 tlie interaction-induced localization wins. This is 
analogous to the case of the surface of a 3D topological insulator discussed above. 
The edge of a 2D topological insulator is protected from the full localization. ^^^ 
This means that the topological distinction between the two insulating phases (or- 
dinary and QSH insulator) is not destroyed by the interaction, whereas the super- 
metallic phase separating them disappears. Therefore, the transition between two 
insulators occurs through an interaction-induced critical state, ^^° see Fig. 1.2 (right 
panel) . 

1.7. Summary 

Despite its half-a-century age, Anderson localization remains a very actively devel- 
oping field. In this article, we have reviewed some of recent theoretical advances 
in the physics of Anderson transitions, with an emphasis on manifestations of crit- 
icality and on the impact of underlying symmetries and topologies. The ongoing 
progress in experimental techniques allows one to explore these concepts in a va- 
riety of materials, including semiconductor structures, disordered superconductors, 
graphene, topological insulators, atomic systems, light and sound propagating in 
random media, etc. 

We are very grateful to a great many of colleagues for fruitful collaboration and 
stimulating discussions over the years of research work in this remarkable field. The 
work was supported by the DFG - Center for Functional Nanostructures, by the 
EUROHORCS/ESF (IVG), and by Rosnauka grant 02.740.11.5072. 
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